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Nondissipative force acting on a vortex in Fermi superfluids and superconductors is calculated 
from the spectral flow of fermion zero modes in the vortex core. The spectral flow effect casts a 
new light on the problem of mutual friction in superfluids. We demonstrate that, for vortices in an 
isotropic system, the reactive mutual friction parameter D'(T) is negative in the coUisionless regime 
. . . and it is positive in the hydrodynamic regime. This agrees with 3He-B measurements [1] ( T.D.C. 

' Bevan, et al, Phys. Rev. Lett., 74, 750 (1995)), which thus give the experimental verification of 

the spectral flow effect in the vortex dynamics. General expression for all 3 nondissipative forces 
(Magnus, lordanskii and spectral flow) is presented in the whole temperature range. 

p PACS numbers: 67.57.Fg, 74.60.Ge, 11.27.+d 
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Spectral flow of fermions is the topological phenomenon related to anomalies in the relativistic quantum field 
theory. The inhomogeneous vacuum induces a spectral flow, which carries particles from negative energy states of the 
vacuum to positive energy states. As a result, the corresponding conserved quantity (viz, charge) is transferred from 
' ^ ' a coherent vacuum motion into incoherent fermionic degrees of freedom, which is visualized as creation of charge from 
■ the vacuum. It was shown in |^,^ that the same phenomenon takes place in superfluid '^He-A, where the gap in the 
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quasiparticle spectrum has topologically stable point nodes, resulting in a spectral flow of Bogoliubov quasiparticles 



0^ through the nodes when the condensate evolves in time. This leads to transfer of a linear momentum from the 
coherent superfluid motion into the normal component of the liquid, i.e., to a force between the superfluid and normal 
0^ ■ components. 



A mutual friction force of the same type arises for a Fermi superfluid or superconductor even without gap nodes, 
if quantized vortices are present 0|. The singularity at the vortex axis plays the same role as gap nodes in A-phase 
^ of ■^He, while the evolution of the condensate, which leads to the spectral flow, is realized by the motion of the vortex 
^ \ with respect to the heat bath. A microscopic analysis of the spectral flow in vortices was also presented in but 
^ ] the hydrodynamic and coUisionless regimes were not distinguished there which led to the contradictory conclusion. 
' Here we consider the influence of the spectral flow on a low-dissipative vortex dynamics in both regimes for isotropic 
' Fermi systems such as s-wave superconductors and '^He-B. For superconductors, we consider a simple case of Galilean 
^ invariant system, i.e., we ignore energy-band effects in metals and assume that the mean free path is longer than the 



coherence length (clean case). We compare the theoretical mutual friction force in '^He with the experimental results 
of Ref . . The observed change of sign of the parameter D' following the change of regimes illustrates the essential 
contribution of the spectral flow effect into the vortex dynamics. 

Nondissipative forces in Bose superfluids. There is a big difference between dynamics of quantized vorticity in 
Bose and Fermi superfluids. In Bose liquids the nondissipative force acting on the vortex moving with respect to 
superfluid and normal components of the liquid is simply a sum of the Magnus forces from the superfluid and normal 
components: 



nd 



KZ X Ps{T){-Vl - V^) +KZX p^{T){vl - v„). (1) 



Here is the vortex velocity, and v„ are the superfluid and normal velocities; kz is the circulation vector. For 
Bose superfluids n — 2TrNh/m, where N is the vortex winding number, and m is the bare mass of a boson. The 
densities Ps{T) and Pn{T) of the normal and superfluid components, respectively, can be anisotropic tensors. 
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Eq. (|l|) can be rearranged to comply with the Landau picture of a superfluid, where the motion of a superfluid 
vacuum of the total density p with the superfluid velocity Vs is supplemented with the dynamics of elementary 
excitations moving with the velocity v„ — . The force becomes 

Fnd = KZX p{vL - Vs) + KZ X Pn{T){^s - V„) 

F]V[agnus Fjordanskii (2) 



Now the Magnus force is deflned as if the vortex moves with respect to the superfluid vacuum. The lordanskii force 
results from elementary excitations outside the vortex core: the vortex line provides them with an Aharonov-Bohm 
potential. 

The lordanskii force in Eq. (||) is a sum of elementary forces which act on individual particles, (9tp = (V x Vs) x p. 
Here p is the quasiparticle momentum and the vorticity V x = K.i62 (r) is concentrated in the vortex core " tube" . 
One has 



X. f^tP = - J nixp f[Ep + p • (vs - v„)] 



(27r3) 

f 

= KZ X p„(T)(v„ - Vs). (3) 

Here f{E) is the distribution function of elementary excitations, Doppler shifted due to the counterflow v„ — v,. The 
lordanskii force has the same origin as the Aharonov-Bohm effect for spinning cosmic strings . 

Adding a dissipative friction force to Eq. (|l|), one obtains a conventional expression for the balance of forces acting 
on the vortex [|o| 

Kps{-Vs - Vl) X z + D(v„ - Vl) - D'{Vn - Vl) X z = 0. (4) 

Here the term with D is associated with the friction while the term with D' is due to the reactive force Fnd- Comparison 
with Eq.(|) gives the reactive parameter 

D' = -Kp„(r). (5) 

Our consideration is valid when dissipation can be neglected compared to the reactive terms, D <g D'. It is only 
under this condition the nondissipative forces do not depend on details of the core structure and of the quasiparticle 
kinetics, and are defined by general characteristics of the vortex and the superfluid. Eq.(Q) coincides with microscopic 
calculations (see review pT|]). 

Nondissipative forces in Fermi superfluids. Now we consider Fermi systems such as superfluid phases of '^He and 
superconductors. For the latter we assume a large magnetic-fleld penetration depth compared to the intervortex 
distance: in this case the magnetic field can be treated as homogeneous. We also assume that the magnetic induction 
B is much smaller than the upper critical field Hc2, so that vortex cores do not overlap. 

In Fermi superfluids, there is a new contribution to Fnd in addition to Eq. (j^) (now k = TrNTi/m and m is the bare 
mass of a fermion): 

Fnd — F]V[agnus ~t~ -^lordanskii -^sp.fiow 7 (6) 

This contribution is related to the chiral fermion zero modes, which are absent in Bose superfluids. 

The spectrum of single-fermionic excitations in a vortex, En(j)z,Q), depends on the momentum projection pz on 
the vortex axis; the orbital quantum number Q, integer or half of odd integer, which corresponds to the generalized 
angular momentum conserved in an axisymmetric vortex; and n denotes radial and spin quantum numbers. The 
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inter level distance of bound states dEn/dQ = ujn is small compared to the gap amplitude A(T) of fermions in bulk: 
LUn (T) /Ep. This spectrum has anomalous (chiral) branches Eq (pz , Q) whose number N^m is related to the vortex 
winding number N^m = '2N according to the index theorem As a function of Q, each anomalous branch crosses 
zero of energy an odd number of times and runs through both discrete and contunuous spectrum from Eq — ~oo 
to Eq — +00. Any other branch either does not cross zero of energy at all or crosses it an even number of times. 
For low-energy bound states, the spectrum of the chiral branch is linear in Q. For the most symmetric vortices, for 
example, Q 

EoiPz,Q) = Qu;oiPz). (7) 

Due to an odd number of crossings of zero, the chiral branch can result in a momentum exchange between fermions 
localized in the vortex core and fermions in the heat bath. The corresponding contribution to the force thus depends 
on v„ — v/,. In contrast to the two other forces, this contribution does depend on the core structure and on the 
quasiparticle kinetics: It is determined by the parameter loqt, where r is the lifetime of fermions. Two limiting cases 
are important, the collisionlcss and hydrodynamic limits: they correspond to loot ^ 1 and ujqt <^ 1, respectively. In 
both cases dissipation is small [[T^-^Tst and the results should thus not be sensitive to the core structure of the vortex. 

Hydrodynamic limit. In the hydrodynamic limit luqt <^ 1 the interlevel spacing is smaller than the level width 1/r 
and the spectrum Eq(pz, Q) can be considered as a function of the continuous parameter Q; it crosses zero at some 
value of Q [at Q = for the spectrum of Eq. When the vortex moves with respect to the heat bath (normal 

component), the velocity difference v„ — induces a flow of quasiparticles from negative levels to positive levels of 
the spectrum Eo{pz, Q). Since the angular momentum evolves as Q ^ Q + (r(i) x p) • z = Q + t((vi — v„) x p) • z, 
the number of levels crossing zero energy per unit time is dtQ = (v^ — v„) • (p x z). Each level bears the linear 
momentum p, therefore, the total flux of the linear momentum from the vortex to the heat bath is 

n.Q 



PF 



dpz 
2ti 



ttTvI^zx (vi-v„). (8) 



Here we used that only zero modes contribute the sum — J^n f i-^n) I ^Q) — Tlmifi-^n ~ ^°°)^ fi^n = 00)) = 2iV 
as Q together with Eo{Q) run from —00 to +00. Thus the spectral flow contribution to F„d in the hydrodynamic 
limit is 

Fsp.flow = KZX Co(v„ - Vi), UJqT <C 1, (9) 

where the anomaly parameter Cq is expressed in terms of the Fermi momentum pp and coincides with the mass 
density of the normal state: 

Co^mp^p/Sn^ . (10) 

The spectral flow of the fermion zero modes within the vortex core is a realization of the Callan-Harvey mechanism 
of the anomaly cancellation in the relativistic quantum field theories [^,^ . 

Collisionless limit in uncharged Fermi superfluids. In the collisionlcss limit loqt ^ 1 the level width is much smaller 
than the interlevel distance and the spectral flow along the discrete levels is suppressed. The spectral flow becomes 
possible only for Eq {p^ , Q) above the gap in the region of the continuous spectrum. Since only the levels in the range 
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A(T) < Eo{pz, Q) < oo are "fluid" the spectral flow takes place only at nonzero T due to the thermal tail of Fermi 
function /(£;„) in Eq.(|). It is reduced by the factor [/(-A(r)) - f{-oo)] - [/(A(T)) - /(oo)] = 1 - tanh[A(r)/2r] 
as compared to the hydrodynamic limit (here we assume that the gap is isotropic). 

Thus for an isotropic pair-correlated system, such as '^He-B, which has an isotropic gap A(T), 



sp.fiow 



— kCo 



1 — tanh 



A(r) 



2T 



z X (v„ - vl), luqt > 1. 



(11) 



Spectral flow in superconductors. For a charged Fermi-system the levels above the gap are also discrete due to 
quantization in the magnetic field, the interlevel distance being the cyclotron frequency lUc- Note that Wc ^ t<^o due 
to the condition B <C Hc2- The levels can be considered as a continuum only imder the condition uJcT 1. The 
results Eqs.(^ and ( pT| ) are thus reproduced for superconductors in the limit lUcT <C 1. In the extreme collisionless 
limit, when lUc-t ^ 1, the spectrum is discrete everywhere and the spectral flow is completely suppressed. Hence, 



^sp.fiow = hCqz X (v„ - Vi), cvqt < 1, ujcT < 1; 

vl), 



A(T)) 

Fsp.flow = kCo (1 - tanh ^„ ) z x (v. 



2T 

wqt > 1, ujcT -c 1; 

Fsp.fiow = 0, UqT > 1, LOcT > 1. 



(12) 

(13) 
(14) 



Parameter D' for ^He-B. Let us consider first the effect of the spectral flow on the mutual friction parameter D' 
in ■^He-B. The dissipative force, proportional to D, is small in both limits, ujqt <^ 1 and luqt ^ 1 [Q. Comparing 
Eq.(^ with Eqs.(||,^,^, we find that the reactive parameter in the hydrodynamic limit is 



Di 



hydro 



= k[Cq - p„(T)], ujoT < 1- 



(15) 



Here it is used that the normal density p„ in '^He-B is isotropic. In the collisionless limit we find from Eq. (|ll| ) 



Co(l-tanh^) -p„(T) 



, LOOT > 1- 



(16) 



Eq.(|l5|) was obtained in Q in the limit T ^ and was extended to nonzero T in HJT^. In the weak coupling 
approximation the anomaly parameter Cq ~ p, since {p — Co) ~ pA'^{T)/Ep ^ p. 
If one neglects the difference between Co and p, 



PsiT) - ptanh 



D'hydToi^PsiT), ujoT < 1, 
A(T)- 



2T 



UJOT > 1- 



(17) 
(18) 



Eq. (^7|) in the limit T <^Tc was derived for superconductors in ||l^ and for ^He in it was extended for higher 
temperatures in ||l^. Eq. (|l^) was first obtained for superconductors in ||l8|,|9| by calculating the force produced by 
excitations scattered by moving vortex. The same result has been reproduced in [l5|] using the kinetic equation. The 
classical results of Bardeen and Stephen and of Nozieres and Vinen |2|] can also be compared with Eqs. (|l^) and 
( p^ ) : The former corresponds to the hydrodynamic limit D' = nps . The result of |2^ , namely that the only force left 
at T = is the Magnus force, remains valid in a collisionless limit, when both the lordanskii force and the spectral 



flow disappear. Its generalization which extrapolates between the two limits |22 does not contain the spectral flow 
term 1 - tanh[A(T)/2r] and thus is vahd only at T = 0. 

Applying the Eqs.( |l7| , |l8| ) to "^He-B one flnds that these two regimes take place in the regions T <^Tc and Tc—T ^ Tc. 
For T close to T^,, the hydrodynamic limit is always realized since Wq '~ A'^(T)/ Ep — > 0, and the reactive parameter is 
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positive D' k, Kps{T). At low T ^ A(0) the collisionless limit is always reached, hence the parameter D' is negative 
according to Eq.(^8|). Indeed, the ratio of D' values for these two regimes is 

^coii/^hydro = 1 - ^ tanh ^ < . (19) 

This dependence is shown in Fig.l. The sign reversal of D' is clearly seen in the oscillating membrane experiment in 
3He-B 0. 

Using the results of Ref . , we can write a phcnomcnological expression which interpolates between the hydro- 
dynamic and collisionless limits: 



I?' « K 



^0 - Pn{T) - ^0 ^^^^^ ^ • (21) 

1 + WgT^ Zl 



The temperature dependences of D{T) and D'{T) can be visualized using a simple model expression for the lifetime 
T ^ (Ep/T^) exp[A(T)/r], so that the parameter luqt is approximated by 

c^or = a[A(r)/T]2exp[A(r)/T], (22) 

where a is a fitting constant. The plot in Fig. 2 qualitatively agrees with the experimental temperature dependence 
found in 0]. We see that D'/nps approaches its negative coUisionless-limit asymptote of Eq.(p^ at low T and 
tends to its positive hydrodynamic-limit asymptote D' / npa = 1 at high T . The observed sign reversal of D' is the 
experimental illustration of the spectral flow force in the vortex dynamics. Eq. (|2^) suggests that for T ^ 0, one 
always has D' « ~npn since Pn/ P ^ 1 — tanh[A(0)/2r], and {uj^t)^"^ decreases faster than p„. 

Superconductors. Eq. ( pl| ) can also be used for investigation of the sign of the Hall effect in superconductors. The 
Hall conductivity in the limit ^ 1 is 



aH^c^{^^.Ps-D')/'^'aB 

ec 



mB 



Co + Co,^tanh^(^) 



1 + luIt^ 2T 



(23) 



where $o — ttc/c is the magnetic flux quantum. If Co > p the sign reversal of an can occur in the hydrodynamic 
regime as proposed in Ref. Q . 

Conclusion. We found a spectral flow contribution to the nondissipative force acting on a moving vortex in 
hydrodynamic and collisionless limits. Due to the spectral flow the reactive mutual friction parameter D' has different 
signs in these two limits. This sign reversal was observed in the ^He-B experiments on the transition between the two 
regimes, realized at different temperatures. This conflrms the essential role of the spectral flow effect in the vortex 
dynamics of Fermi superfluids and supercoductors. 
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FIG. 1. Theoretical dependence of the ratio D'^oiil ^[^y^,.^ from Eq. (^|) for an isotropic superfluid. This ratio is always 
negative. The gap function A(r) and the superfluid density PeiT) are taken in the weak-coupling BCS approximation. F\ is 
the Fermi-liquid parameter which enters pa in a Fermi superfluid. For ^He, the values of Fi_ correspond to pressures p = Q, 
p = 12, and p = 24 bar, respectively. 



FIG. 2. Qualitative temperature dependences of D/hips and D' /ups from Eq. ( |2l[ ) in a simple model of Eq.(|22[). The fitting 

parameter a = 0.02. 





6 



